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We apply the Dijkstra algorithm to generate optimal paths between two given sites on a lattice 
representing a disordered energy landscape. We study the geometrical and energetic scaling prop- 
5^ ' erties of the optimal path where the energies are taken from a uniform distribution. Our numerical 

Oh. results for both two and three dimensions suggest that the optimal path for random uniformly dis- 

tributed energies is in the same universality class as the directed polymers. We present physical 
realizations of polymers in disordered energy landscape for which this result is relevant. 
Recently, there have been much interest in the prob- 
lem of finding the optimal path in a disordered energy 
r TT 1, landscape. The optimal path can be defined as follows. 
'■^ \ Consider a d-dimensional lattice, where each bond is as- 
■ signed with a random energy value taken from a given 
distribution. The optimal path between two sites is de- 
I ■ fined as the path on which the sum of the energies is 
| minimal. This problem is of relevance to various fields, 
such as spin glasses Q , protein folding ||] , paper rupture 
H , and traveling salesman problem 0] . Though much ef- 
d • fort has been devoted to study this problem, the general 
solution is still lacking. There exist two approaches de- 
ll ■ veloped recently to study this problem. Cieplak et al. || 
applied the max-flow algorithem for a two-dimentional 
q ■ energy landscape. Another approach is to restrict the 
O \ path to be directed, that is, the path can not turn back- 
wards. This approach is the directed polymer problem 
I ' which has been extensively studied in the past years, see 
e.g., §§■ 

' In this manuscript, we adapt the Dijkstra algorithm 
from graph theory || for generating the optimal path on 
a lattice with randomly distributed positive energies as- 
. signed to the bonds. This algorithm enables us to gener- 
Q\ ' ate the optimal path between any two sites on the lattice, 
ON . not restricted to directed paths. We study the geometri- 
cal and energetic properties of the optimal paths in d = 2 
and 3 dimensions in random uniform distribution of en- 
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ergies. We calculate the scaling exponents for the width 
and the energy fluctuations of the optimal path. We find 
that for both d = 2 and d = 3 the exponents are very 
close to those of directed polymers suggesting that the 
' non-directed optimal path (NDOP) is in the same uni- 
J> . versality class as directed polymer (DP). Our results are 
in agreement with those found by Cieplak et al. fa] for 
the two dimentional case. This result indicates that in 
the case of uniformly distributed energies the NDOPs are 
self affine and overhangs do not play an important role 
in the geometry of NDOPs. 

Our results are relevant, for example, in the following 
polymer realizations: (i) Consider a d-dimension energy 
landscape in which there is a spherical regime of ran- 
domly distributed high energies while outside this sphere 
the energies are zero or have very low values. Consider 
as well a polymer of length N that one of its ends is 
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attached to the center of the sphere while the other is 
free. The radius of the sphere is r <C N (see Fig la). 
The section of the polymer inside the sphere will reach 
the lowest energy path which is the optimal path stud- 
ied here i.e., with a self afhne structure, (ii) Consider 
a polymer in a d-dimensional energy landscape which is 
divided into alternating strips of disordered low and high 
energies (see Fig lb). In the strips of high energies the 
polymer is expected to behave like the optimal path. 

The Dijkstra algorithm enables one to find the opti- 
mal path from a given source site to each site on a d- 
dimensional lattice. During the execution of the algo- 
rithm, each site on the lattice belongs to one of three 
sets (see Fig. ||): 

• The first set includes sites for which their optimal 
path to the source site have been already found. 

• The second set includes sites that are relaxed at 
least once but their optimal path to the source has 
not been determined yet. This set is the perimeter 
of the first set. 

• The third set includes all sites on the lattice which 
have not been visited yet. 

The algorithm itself consists of two parts, (i) initializa- 
tion and (ii) the main loop. The main loop, in its turn, 
is composed from (a) the search and (b) the relaxation 
processes. 

In the initialization part we prepare the lattice in the 
following way. Each bond is assigned with a random 
energy value taken from a given distribution. Each site 
is assigned with an energy value of infinity. We pick up 
a certain source site and assign it with the energy value 
of zero and insert it into the second set. 

After that we enter the main loop. We perform the 
search among the sites from the second set and find the 
one with the minimal energy value. Then we add it to the 
first set and proceed to the relaxation process. This site 
is called the added site. The relaxation process deals with 
sites neighboring to the added site that do not belong to 
the first set. 

In the relaxation process we compare two values, the 
energy value of the neighboring site and the sum of two 
values: the energy value of the added site and the energy 
value of the bond between these sites. If the value of the 
sum is smaller, then we a) assign it to the neighboring 
site, b) connect the neighboring and the added site by 
the path (thick bond in the Fig. ^) , c) if the neighboring 
site belongs to the second set we break its previous con- 
nection to another site (thick bond) , d) if the neighboring 
site does not belong to the second set we insert it into 
the second set. The first four steps are demonstrated in 
Pig. | 

Normally, the main loop stops when the second set is 
empty, however one might wish to break the loop earlier, 
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e.g., at the moment when the first set reaches the edge 
of the lattice in order to avoid boundary effects. 

Each site that belongs to the first set is connected to 
the source by a permanent path (thick bonds) that does 
not change during the execution of the algorithm, so, if 
we stop the algorithm at any given time the first set will 
still be valid. 

We simulate both DPs and NDOPs on a square lattice 
in the following way. Let x, y be the horizontal and ver- 
tical axes. We choose the origin to be the source site and 
study the optimal paths connecting it with all the sites 
on the line between [0, t] and [t, 0] for different values of 
t. The generalization to three dimensions is straightfor- 
ward. The random energies assigned to bonds are taken 
from a uniform distribution. We find that our results are 
independent of the distribution interval. 

In Fig. H we compare a configuration of DP and NDOP 
on the same disordered energy landscape. It is seen that 
in the NDOP only very few overhangs exist. To test the 
effect of the overhangs we calculate the mean end-to-end 
distance R of the global optimal path (thick line in Fig. ||) 
as a function of its length I. The global optimal path is 
the minimal energy path among all the paths with the 
same value of t. Our numerical results clearly indicate 
the asymptotic relation I ~ R showing that the NDOPs 
are self-affine 0. We should compare this result to the 
strong disorder limit where the paths can be regarded as 
self-similar fractals, with I ~ i? dopt , where d op t — 1.22 in 
d = 2 and d opt ~ 1.42 in d = 3 0,0. 

In order to compare NDOP to DP we study several 
properties, such as the roughness exponent £, the en- 
ergy fluctuation exponent £ for two and three dimen- 
sions, as well as the distribution of the endpoints of DP 
and NDOP. The above exponents are defined by the rela- 
tions W = (li 2 ) 1 / 2 ~ t* and AE = ((E - (E)) 2 } 1 / 2 ~ t*. 
Here, h is the transverse fluctuation of the global optimal 
path which is distance between its endpoint and the line 
x = y; E is the energy of the global optimal path which is 
the sum of all bond energies along the path. The average 
is taken over different realizations of randomness. Fig. |] 
shows the dependence of width W and energy fluctuation 
AE of DP and NDOP on t in two and three dimensions. 
The points are the data for both DP and NDOP and the 
dashed lines represent the exponents of the DP. Our re- 
sults indicate that the exponents for the NDOP are very 
close to those of DP (see also Table |) . 

Our results may be related to the recent findings ]T^1 
that the roughness exponent of the minimal energy of 
domain wall in random Ising model and fracture interface 
are the same H] in d = 2. In these cases, similar to our 
case, although overhangs may occur they do not play an 
important role. 

In summary, our results suggest that the optimal path 
in the case of uniformly distributed energies, for any en- 
ergy interval, is in a different universality class from the 
strong disorder limit but in the same universality class 
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as directed polymers. This result is relevant to several 
questions regarding the equilibrium state of polymers in 
different realizations of disordered energy landscape. 
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TABLE I. Width and energy fluctuation exponents of DP and NDOP in 2 and 3 dimensions. The exponents were derived 
from the slopes of the corresponding data points shown in Fig. ^ The error bar were estimated from taking five ensembles of 
10 4 configurations each for d — 3 and 500 configurations each for d — 3. 
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NDOP 




0.66 ± 0.02 


0.67 ±0.02 


0.60 ± 0.05 


0.63 ± 0.05 
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0.32 ± 0.02 


0.32 ±0.02 


0.19 ±0.07 


0.18 ±0.07 
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FIG. 1. Schematic illustrations of sections of a polymer having the structure of optimal path. (a) within the circle energies 
are distributed and outside the circle zero or very low energies are distributed, (b) alternating strips of high and low distributed 
energies. 
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FIG. 2. Illustration of first four steps of the Dijkstra algorithm applied to a square lattice. Numbers along the bonds 
represent the random energy assigned to them. Numbers inside circles represent energies of sites i.e., the total energy of the 
path connecting this site to the source. Empty sites possess infinite energy and belong to the third set, thick circles belong to 
the first set, all other circles belong to the second set. Note e.g., that the site marked 6 in (c) got relaxed one more time and 
become 5 in (d) . During this second relaxation we break its previous connection to the site with energy 2 and connect it to the 
site with energy 3. At each time step we identify the optimal path from each site in the first set to the source by going along 
the thick bonds. 



G 



FIG. 3. The sets of all directed (the upper one) and non-directed (the lower one) optimal paths with t = 300 obtained for 
the same realization of quenched randomness in the lattice. The global optimal path which is the minimal energy path among 
all the paths with the same t is shown by a thick line. In this particular case the directed and non-directed global optimal 
paths do not overlap. In other cases they might overlap significantly though the rest of tree looks somewhat different. 
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Fig. 4a 
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FIG. 4. The width and the energy fluctuation as a function of t on a double logarithmic plot in (a) two and (b) three 
dimensions. Circles are used for directed polymers and squares for non-directed optimal paths. For (a) 10 systems of linear 
size up to t = 300 and for (b) 5000 systems of linear size up to t = 75 are used. The dashed lines are given as a guide to the 
eye and have slopes equal to the exponents known for DP: £ = 2/3, £ = 1/3 for d — 2 and £ ~ 0.62, £ ~ 0.24 for d = 3 |J. For 
all cases we used a uniform distribution of energies between Ei — 1 and E2 = 1000. We also tested other energy intervals and 
found the same results. 
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